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SmiMARY 


A local inviscid-viscous interaction technique has been developed for the 
analysis of low-speed airfoil leading-edge transitional separation bubbles. In 
this analysis an inverse boundary-layer finite-difference analysis is solved 
iteratively with a Cauchy integral representation of the inviscid flow which is 
assumed to be a linear perturbation to a known global viscous airfoil analysis. 
Favorable comparisons with data indicate the overall validity of the present 
localized interaction approach. In addition numerical tests were performed to 
test the sensitivity of the computed results to the mesh size, limits on the 
Cauchy integral, and the location of the transition region. 


2 



INTRODUCTION 


The classification of airfoil stall has been well known since the investigation 
of Jones (Ref. 1) in the early 1930 *s showed that there are three -basic types of 
stall at subsonic speeds induced by trailing-edge separation and two types of 
leading-edge separation. Later the detailed experiments of McCullough and Gault 
(Ref. 2) revealed that trailing edge stall is associated with the separation of 
a turbulent boundary layer which moves forward with increasing incidence, 
leading-edge stall results from an abrupt separation of the laminar boundary layer 
near the leading edge without subsequent reattachraent , and thin airfoil stall in 
wMch the laminar boundary layer separates near the leading edge followed downstream 
by turbulent reattachment which progresses rearward with increasing incidence. 

Further discussion on .the classification of airfoil stall is presented in the 
excellent review article by Tani (Ref. 3). 

In addition there has been recent work by van de Berg (Ref. A) which indicates 
that leading-edge stall may not be the result of bursting of the laminar separation 
bubble, but rather it is due to separation of the turbulent boundary layer which 
develops downstream of this bubble. In this case it is still necessary, as pointed 
out by van den Berg, to be able to analyze the flow in the laminar separation 
bubble since it forms the initial conditions of the turbulent boundary layer 
downstream. 

The focus of the present investigation is the development of a prediction 
technique for closed airfoil leading edge separation bubbles. It was pointed out 
by Tani (Ref. 3) that airfoils at moderate incidence angles, prior to either 
leading-edge stall or thin airfoil stall, experience local separation bubbles just 
downstream of the peak suction (minimum pressure) region. Figure 1 shows a 
schematic diagram of an airfoil leading-edge bubble which occurs if the Reynolds 
number is sufficiently low so that the boundary layer remains laminar up to the 
minimum pressure point. Downstream of this point separation occurs almost 
immediately since laminar boundary layers, in contrast with turbulent flows, are 
extremely sensitive to adverse pressure gradients. A separation bubble forms in 
which a recirculating streamline pattern is bounded by a shear layer. Since shear 
layer flows tend to be highly unstable to flow disturbances, transition from laminar 
to turbulent flow generally occurs in this shear layer. Further downstream the 
turbulent mixing between the shear layer flow with the lower dead air region results 
in entrainment of higher energy air which energizes the flow near the surface thereby 
resulting in flow reattachment with subsequent turbulent boundary-layer flow down- 
stream, As shown in Fig. 1 the initial portion of the separation bubble is 
characterized by a pressure plateau followed by a pressure recovery region after 
the transition process is initiated, but prior to flow reattachraent. An increase 
in incidence causes the bubble to move forward and contract in strearawise extent 
until the flow no longer reattaches. At this incidence angle bubble bursting has 
occurred, thereby resulting in leading-edge stall which is accompanied by an abrupt 
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loss of lift since the suction peak has now collapsed with the resultant pressure 
distribution redistributed in a flattened form over the airfoil chord. 

There have been numerous experimental investigations conducted such as the 
work of Bursnall and Loftin (Ref. 5), Gault (Ref. 6), Caster (Ref. 7), and more 
recently that of Mueller and Batill (Ref. 8) to provide information on the flow in 
the neighborhood of laminar transitional separation bubbles. Owen and Klanfer 
(Ref. 9) deduced from experimental measurements that bubble bursting occurs if 
the momentum thickness Reynolds number at the laminar separation point is less 
than 125; similarly, Crabtree (Ref. 10) proposed a criterion based on the pressure 
rise over the bubble. Later, Horton (Ref. 11) developed a semi-empirical theory 
based on the experimental measurements of Caster (Ref. 7) for the growth and 
bursting of laminar separation bubbles. At the present time airfoil analysis 
codes such as the NASA-Lockheed multielement airfoil code (Ref. 12), and the 
GRUMFOIL code (Ref. 13) use simple criterions such as these to detect the 
occurrence of laminar separation bubbles and whether or not bursting occurs. 

In these analyses if laminar separation is detected without bursting then the 
flow is assumed to immediately undergo transition to turbulent flow with a jump 
discontinuity in the boundary-layer parameters such as shape factor and skin 
friction. 

With the recent theoretical developments in boundary-layer and viscous-inviscid 
interaction theory for separated flow, there have been several analytical investigations 
(Refs. 14-18) conducted to provide a more detailed description of the flow process in 
a laminar separation bubble. These studies have been based on both finite-difference 
and integral techniques for describing the boundary-layer flow which are solved 
iteratively with a Cauchy integral representation of the inviscid flow. A particularly 
encouraging result was obtained by Briley and McDonald (Ref. 14) in which they showed 
that solutions obtained with viscous-inviscid theory were essentially the same as 
solutions obtained with the more computationally expensive Navier-Stokes equations. 

The overall objective of the present investigation is the development of an 
analytical procedure for predicting airfoil leading-edge laminar separation 
bubbles. Prior to bursting, the strearawise extent of the separation bubble is 
only a few percent chord; hence, in a manner similar to the previous theoretical 
approaches discussed above, a local viscous-inviscid interaction analysis has been 
developed for the immediate vicinity of the separation bubble. However, in contrast 
with these previous efforts, the present perturbation approach is based on the use 
of a global airfoil viscous analysis which provides the reference pressure and 
displacement thickness as input for the local leading-edge interaction analysis. 

In this manner the global flow field around the airfoil, including the wake, is 
taken into account. In the present study the viscous airfoil analysis of Melnik, 

Chow and Mead (GRUMFOIL code. Ref. 13) has been used. 

Computations are presented for several airfoils for which experimental data 
was available. In addition, comparison is made with the benchmark experimental 
data obtained by Caster (Ref. 7). Overall the computed solutions are in good 
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agreement with the data, although it was observed that the solutions are quite 
sensitive to the transition model as will be demonstrated in the present report. 
Since the major thrust of the present work was to demonstrate the feasibility of 
a local interaction model used in conjunction with a global airfoil analysis, 
a simple streamwise intermittency function has been used to represent the transition 
from laminar to turbulent flow. 
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SYMBOLS 


c airfoil chord 

skin friction coefficient 
f perturbation stream function 

F velocity ratio, u/u 

e 

g total enthalpy ratio, H/H 

e 

H total enthalpy 

L reference length 

M Mach number 

m perturbation mass flow 

n coordinate normal to reference displacement surface 

N coordinate measured normal to reference displacement surface from the body 

surface 

Pr Prandtl number 

Pr^ turbulent Prandtl number 

q perturbation source strength per unit length 

Re reference Reynolds number 

Re^ local momentum thickness Reynolds number 

s, S coordinates along reference displacement surface 

u velocity component parallel to reference displacement surface 

V velocity component normal to reference displacement surface 

V transformed normal velocity in Prandtl transposition theorem 

3 pressure gradient parameter 

V streamwise intermittency function 


6 



6" displacement thickness 

G eddy viscosity coefficient 

n transformed normal coordinate 

y molecular viscosity coefficient 

5 transformed tangential coordinate 

p density 

stream function 
Subscripts 

e edge of boundary layer 

I inviscid 

ref reference solution 

t^ start of transition 

t end of transition 

2 

V viscous 

“ free stream 

1 start of interaction region 

2 end of interaction region 
Superscripts 

’ perturbation quantity 

i global viscous-inviscid iteration counter 
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VISCOUS-INVISCID INTERACTION ANALYSIS 


The approach taken in the present effort to develop a prediction technique for 
airfoil leading-edge separation bubbles is based on a viscous-inviscid interaction 
technique in which the boundary-layer equations are solved iteratively with an 
inviscid analysis through displacement thickness coupling. As mentioned in the 
Introduction, the focus of the current effort is the analysis of closed transitional 
separation bubbles which are known from the experimental studies cited previously 
to occupy only a few percent of the airfoil chord. Since the resultant interaction 
is highly localized, it was decided to treat the leading-edge transitional bubble 
problem as a linear perturbation to a known global airfoil solution. The use of a 
localized interaction analysis, despite the overall approximation inherent in a 
perturbation approach, permits an accurate analysis of the flow field structure in 
this region, in contrast with the extremely difficult problem of trying to resolve 
this small scale phenomena while simultaneously solving the global airfoil flow 
field. However, the present approach, in contrast with previous perturbation 
treatments of this problem, accounts for the influence of the global viscous 
airfoil flow on the local interaction analysis. It is well known that at high 
angles of attack, where transitional bubbles frequently occur, the turbulent 
interaction at the airfoil trailing edge and wake influence the entire airfoil 
pressure including the location of the stagnation point and the level of suction 
pressure, which are the flow conditions upstream of the transitional bubble. In 
the present investigation the GRUMFOIL viscous airfoil analysis, developed by 
Melnik, Chow, and Mead (Ref. 13) , was used to establish the reference pressure and 
displacement thickness distributions which were then input to the local interaction 
analysis which is described in this section. The local inviscid analysis is 
discussed first followed by a brief description of the boundary-layer analysis 
and the iteration procedure used to couple these two analyses. 


Inviscid Analysis 

Since most of the available experimental data for transitional bubbles is for 
low speed flows, the present analysis was restricted to low speed flow for which 
Laplace’s equation is the governing equation for the inviscid velocity potential, $. 
Since this equation is linear, solutions can be superimposed and hence the velocity 
potential is written as 


+ ( 1 ) 

where $ f velocity potential obtained from the global airfoil analysis 

for the flow over the reference displacement surface, shown in Fig. 2, and <()’ is 
the perturbation potential due to the local transitional bubble. In the present 
analysis the disturbance field is represented by a distribution of sources placed 
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along the reference displacement surface for which the solution for <}>* can be 
expressed as 

I r 2 

= q(0 logLls-«> df ( 2 , 

where q is the source strength per unit length and s and n are the coordinates along 
and normal to the reference displacement surface, respectively. Equation (2) is 
only an approximate solution since it assumes that the sources are distributed 
along a straight line, whereas near the airfoil nose the reference displacement 
surface, which closely follows the airfoil shape, is highly curved. However, 
since the interaction length is very short it is assumed that over the interaction 
region the reference displacement surface is a straight line and hence the surface 
curvature has been neglected in the present investigation. 

The perturbation velocity components induced by the deviation of the interaction 
displacement surface, shown in Fig. 2, from the reference displacement surface, are 
given by , 

(3) 

OS 3n 

which can be rewritten in terms of the perturbation potential using Eq. (2): 


u'(s,n) 



(s-e)dC 

(s-C)^+n^ 


(A) 


and 


v'(s,n) 


27t 1 (s-e)2+n‘ 


(5) 


The source strength q is deduced by the imposition of the surface boundary condition 
which requires that the inviscid flow be tangent to the interaction displacement 
surface shown in Fig. 2. This boundary condition is given by 


v'(s.S*)=u(s,S*) 


d(S-Sref) 

ds 


I 

h{s,S*) 


( 6 ) 


where 
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(7) 


u:uref + u' 


and 


h = i + /cn 


( 8 ) 


where < is the curvature of the reference displacement surface which was set to 
zero in the present investigation. The boundary condition given in Eq. (6) can be 
transposed to the reference displacement surface by a Taylor series expansion and 
application of the continuity equation which results in 


v'(s,o) 


d[ue<S*-8re|)] 

ds 


+ o(S*-8%f)^ 


(9) 


where 


Ue(s) = u(s,o)=Uref (s,o) + u'(s,o) 


( 10 ) 


The source strength is deduced by taking the limit of Eq. (5) as n -► 0 and the use 
of Eq. (9) which leads to 


q(s)=2v'(s,o) 


ds 


( 11 ) 


The final result needed for the present viscous-inviscid interaction calculations is 
the Cauchy integral for u'(s,o) which is obtained by the substitution of Eq. (11) into 
Eq. (4) which gives 


u'(s,o) 


i s-e 


( 12 ) 


with the total velocity u (s) , which is tangent to the interacted displacement surface, 
deduced from Eq. (10). 

The Cauchy integral given in Eq. (12) is evaluated from s^ to s^ which is the 
region of strong interaction that results from the presence of the local transitional 
separation bubble. In the present calculations the reference displacement surface is 
computed from the GRUMFOIL analysis in which instantaneous transition from laminar 
to turbulent flow is assumed to occur at the predicted laminar separation point. 
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As a result the interaction displacement surface merges smoothly with the reference 
displacement surface upstream of the interaction as shown in Fig. 2. Hence the lower 
limit, s^, on the Cauchy integral is placed sufficiently far upstream of the interaction 
region where the source strength is zero. The placement of the downstream limit is 
somewhat more complicated since the interacted displacement surface does not merge 
into the reference displacement surface downstream of the interaction region. This 
result is expected since the presence of the transitional separation bubble results 
in a thicker downstream turbulent boundary layer than that which is obtained with 
the GRUMFOIL analysis which eliminates the transitional separation bubble by placing 
transition at the laminar separation point. However it was found in the present 
calculations that the source strength is very close to zero downstream of the bubble 
since the interacted and reference displacement thicknesses are nearly parallel to 
each other. Numerical tests were performed which confirmed the insensitivity of 
the local interaction analysis to the values of the limits s^ and s^ provided they 
were placed sufficiently far enough upstream and downstream, respectively, of the 
interaction region. It is noted that the present perturbation procedure, which is 
based on the use of a reference displacement surface, considerably simplifies 
the choice of these limits over that used by other investigators in which the 
Cauchy integral was used to predict the disturbance field due to the total 
interaction displacement surface. In the latter case the reference solution is 
the inviscid flow over the airfoil. 

The Cauchy integral given in Eq. (12) is evaluated numerically using a second- 
order scheme developed by Napolitano, et. al (Ref. 19). In this procedure the 
integration region, s < s < s^, is subdivided into I ^ grid points with the complete 
integral written as the sum of integrals over these su§intervals to yield the 
perturbation at the ith grid point: 


^max-i 



/ 

Vi 


-^(UeA8*)de 


(13) 


where 


AS* = 8*-8,ef 


(14) 


and 


_ Sj+Sj+, 

®j' 2 


(15) 


With 

"^J-1 


the source strength, d/dC (u A6 ) assumed to vary linearly over each subinterval, 
to s., Eq. (13) is analytically integrated to give 
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This numerical procedure permits the use of a nonuniform mesh distribution which 
was used in the present problem to concentrate grid points near the center of the 
interaction region in order to adequately resolve the high gradient phenomena which 
occurs as the flow undergoes transition from laminar to turbulent flow. The same 
surface grid point distribution was used in both the inviscid and boundary-layer 
analyses thereby avoiding interpolation between these two solutions in the 
interaction calculations. 


Viscous Analysis 

The viscous solution technique used in the present investigation is the 
inverse boundary-layer procedure presented by Carter (Ref. 20) for the analysis of 
separated flows. Although the inviscid analysis discussed in the previous section 
has been limited to low speed flow, the boundary-layer analysis used in the 
present study was adapted from earlier work which was for compressible flow. Thus 
the fully compressible boundary-layer analysis is presented here. 

The nondimens ional boundary-layer equations are written as follows in terms of 
the reference displacment surface coordinate system shown schematically in Fig. 2: 


d p\} 

"IT 



(17) 


du du 


dx 


-A. ^ ^ _ ouv' 
dn dn 


(18) 


pM 


an 

as 


an 


+ />v 


_a 

an 




(19) 


The v-component of velocity and the n-coordinate are scaled in the usual manner by 
VRe where Re is the Reynolds number based on the free stream flow conditions 
and iSe airfoil cfiord. The boundary conditions imposed on these governing equations 
are ; 
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n=-Sjef u=vsO, 


an 

H or 


specified 


( 20 ) 


n— 00 u — Ue, H— He 

The Reynolds stresses are related to the mean flow by 


( 21 ) 


n-r -7-T JJlJiL 

U V =€y • -P v'h' = Prt an 


( 22 ) 


where e is the eddy viscosity coefficient for which the Cebeci-Smith (Ref. 21) two- 
layer model for fully turbulent flows was used. The transition from laminar to 
turbulent flow was modeled by the use of a streamwise intermittency function y(s) 
which varies from 0 to 1 over a specified region. The treatment of the transition 
region will be discussed in more detail in the next section. 

Werle and Verdon (Ref. 22) showed that it is convenient to transform the boundary- 
layer equations with the Prandtl transposition theorem which is given by 


S*s, N?n+ 8 ref » 


V=v + u 


ds 


(23) 


where N is a transformed normal coordinate measured from the airfoil surface 
perpendicular to the reference displacement surface. In Ref. 22 it is shown that the 
form of the boundary-layer equations is unchanged by this transformation and the 
same surface boundary conditions (Eq. (20)) are imposed at N = 0. 

The development of the inverse formulation begins by transforming the equations, 
expressed in primitive variables, by the following transformation of the independent 
variables 


X 

0 



(24) 


which is quite similar to the Levy-Lees transformation. It is helpful to scale the 
normal coordinate by the displacement thickness in strongly interacting flows since 
this step insures that the boundary- layer thickness is approximately constant in the 
transformed coordinate. The continuity equation is eliminated by introducing the 
stream function 
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The value of the stream function at the boundary-layer edge is written in terms of 
the displacement thickness 


'I' Ug( N- 8*) 

where 



(26) 


(27) 


Then with the definitions 


m 



o 


the edge value of the stream function can be written as 


(28) 


m {77 - 1 + h) as 77—^00 (29) 


A perturbation stream function is defined as 


7 = ^ [+-Fm (-)-l+h)] with F = ( 30 ) 

such that f 0 as n ~ for a prescribed m. Note that in the transformed inverse 
formulation the perturbation mass flow m is prescribed and not just the displacement 
thickness. 

Transformation of the compressible boundary-layer equations with Eq. (24) and 
introduction of the perturbation stream function defined in Eq. (30) results in the 
following set of governing equations: 



(31) 


AT - m 


(I - 17 - h) 


dr) 




m^)3{g-F2) 



(32) 


m^F -m^ [y^ T+ mF(,, - I + h)]|^ = 

±±\J^^S.fL\^Q] ±L{. npill 

pr dr, t^Pr^ldr)!^ I + id. m 2 dr, [M' Prf^dr)] 


(33) 


where 

a JL R -L.^ i It. 

9'He ^"Med^ PtH-t (34) 

Equations (31) - (33) are solved for F, g, f and 3 for a prescribed streamwise 
distribution of m subject to the following boundary conditions: 


17 = 0 


17 — 00 


F = T=o g=g„ or|i 


specified 


F = g — ^ I ond f 0 


(35) 


These equations can also be solved in the direct mode with 3 prescribed and the outer 
boundary condition f = 0 eliminated* In this case if m is set equal to then the 
usual Levy-Lees formulation is deduced with the only difference being that the normal 
component of velocity has been re-expressed in terms of the stream function. In the 
inverse case the unknown pressure gradient parameter is deduced simultaneously with 
the remainder of the solution. The numerical solution of these equations for both 
the direct and inverse mode is an implicit finite-difference technique which is first- 
order accurate in the stream direction and second-order accurate in the normal 
direction. The details of the numerical scheme are presented in Ref. 23. 
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Transition Model 


The prediction of the transition process as a flow changes from laminar to 
turbulent remains one of the most challenging problems in fluid dynamics despite 
the enormous attention given this problem over the past sixty years. The 
prediction of transition is particularly difficult in the airfoil leading edge 
separation bubble since it occurs in a region where the viscous flow is strongly 
interacting with the local inviscid flow. To the author's knowledge, there have 
been only a few limited efforts, such as the work of van Ingen (Ref. 24), which have 
attempted to analyze the instability process which leads to transition from laminar 
to turbulent flow in a separated shear layer bounding a closed recirculation region. 
Since the objective of the present investigation was the development and 
demonstration of a viscous- inviscid interaction method for the analysis of airfoil 
leading edge transitional separation bubbles, a simple forced transition model in 
which the onset and length of transition were prespecified was used to obtain the 
results presented in this report. This model is described in this section along 
with a brief discussion of several attempts that were made in the early part of 
this investigation to use several existing semi-empirical natural transition models. 


The forced transition model used in the present analysis is based on the 
streamwise intermittency distribution which was established by Dhawan and Narasimha 
(Ref, 25) by correlating a number of experimental studies of flows undergoing 
transition. This intermittency function is written as 


r =i-exp 



(36) 


where s is the streamwise position where transition begins and 
^1 

X=sl -si 

1Y=0.75 IY=0.25 


(37) 


The function y> which multiplies the turbulent eddy viscosity as shown in Eq. (22), 
varies form 0 to 1 as the flow undergoes transition from laminar to fully turbulent 
flow. It was found to be convenient to express X in Eq. (36) in terms of the 
location ®t 2 » transition and was assumed to occur at y = .99. 

With this modification, Eq. (36) becomes 



(38) 


In the present calculations, the values of s^ and s^ were deduced from the 
experimental data and input to Eq. (38) to complete tne specification of y(^)* 

The sensitivity of the computed results to the intermittency function is presented 
in the section on Results and Discussion. 
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In the early stages of this investigation, calculations were attempted using 
the semi-empirical methods of Granville (Ref. 26) and Michel (Ref. 27) to predict 
the onset of transition. However, it was found that the use of these natural 
transition models predicted the onset of transition much further downstream than 
that which occurred in the experiments of Caster (Ref. 7) and Gault (Ref. 6) on 
laminar separation bubbles. The main difficulty in the use of these models for 
the Caster and Gault cases is that the Reynolds numbers in these experiments 
(Re = 2 X 10 ) are at the lower Reynolds number range for which these correlations 
were established for which the data was insufficient to establish an accurate 
correlation. In addition, these correlations are probably not applicable to the 
present study of separated flow since they were established for much milder adverse 
pressure gradients than those encountered in the present study. 


Interaction Iteration Procedure 

The present analysis is based on a viscous-inviscid iteration technique which 
was previously presented by Carter (Ref. 28) and is adapted to the present 
investigation as outlined in Fig. 3. This procedure, which has been referred to 
as a semi- inverse technique by LeBalleur (Ref. 29), combines an inverse boundary- 
layer technique with a direct inviscid analysis via the update procedure shown in 
Fig. 3. The key feature of this iteration procedure is the simple update formula 


m 


i+! 




(39) 


which permits an inverse boundary- layer analysis to be directly linked to an 
inviscid analysis which accounts for displacement thickness effects. Note that for 
simplicity the update formula is shown in Fig. 3 with the relaxation factor, o), 
set equal to one. This interaction iteration procedure, shown in Fig. 3, is general 
and has been used in several other studies (Refs. 16, 22, 30 and 34) in which a 
variety of inverse boundary-layer and direct inviscid solution procedures have been 
employed. It was found by Kwon and Fletcher (Ref. 16) that convergence could be 
accelerated by making several inner-loop passes through the Cauchy integral and the 
update formula with the boundary-layer prediction of the edge velocity , u^, frozen 
atits current global iteration value. This technique was used in the present 
calculations with three inner-loop passes and was found to accelerate the global 
convergence rate by a factor of three with a 50% reduction in computer time as 
compared to calculations made without this inner iteration. 
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RESULTS AND DISCUSSION 


The Cauchy integral and inverse boundary-layer analysis which are used to 

iteratively solve for the local interaction region near a transitional separation 

bubble were combined into a computer code referred to as ALESEP (Airfoil Leading 

Edge Separation) • The code is set up to input the displacement thickness and 

pressure distribution for the reference solution about which the local interaction 

analysis is assumed to be a linear perturbation. In this section, the results 

obtained with this analysis are presented for laminar and transitional separated 

flow over a flat plate with a trough, and comparison is made with the experimental 

data of Caster (Ref. 7), and Gault (Ref. 6) for transitional separation bubbles. 

The sensitivity of this interaction analysis to the transition model, Reynolds 

number, grid size, and limits on the Cauchy integral is demonstrated in these 

results. In all of these calculations, 100 grid points were placed nonuniformly 

across the boundary layer with the minimum step size adjacent to the wall prescribed 

so as to provide an adequate resolution of the high gradient turbulent boundary 

layer which occurs downstream of the separation bubble. These calculations were 

initiated with an assumed linear distribution of the total displacement thickness 

and were continued until the difference between ue and ue was less than 

I V 

0.2%. In general about 20 global viscous-inviscid interaction iterations were 
required. In most cases, the calculations were performed with a relaxation factor of 
unity. 


Flat Plate With Trough 

The present analysis was verified by computing the laminar separated flow over 
a flat plate with a trough which was initially analyzed by Carter and Wornom 
(Ref. 32) followed later by Kwon and Pletcher (Ref. 16), Veldman (Ref. 33) and 
Davis and Werle (Ref. 34) . Since the results from all of these analyses are in 
excellent agreement, this separated flow provides an excellent benchmark case by 
which to test a newly-developed low speed viscous-inviscid interaction analysis. 

In this case, the interaction analysis was used to treat the total displacement 
body and hence Ug and were set equal to zero . The interaction region was 

assumed to occur in the region, 1.0 < x/L ^ 4.0, with the trough centered at 
x/L = 2.5 as shown in Fig. 4(a). In this case, the arclength has been replaced in 
the interaction analysis by the Cartesian coordinate x since they are nearly 
equal for this shallow trough. The displacement thickness contribution to the 
Cauchy integral upstream and downstream of the interaction region was treated by 
assuming that in these regions /x. For this case, a streamwise mesh of 121 grid 
points was distributed over the interaction region with a constant step size of 
x/L = 0.025. 

The results obtained with the present analysis for this case were found to be 
in excellent agreement with those obtained previously by Carter and Wornom (Ref. 31). 
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The results for the edge velocity, skin friction, and displacement thickness are 
shown in Fig. 4. In addition, results are presented which demonstrate the effect 
of an imposed transition region, 2.35 < x/L < 2.55, over which the intermittency 
factor was assumed to vary linearly from zero to one. In this case, the displacement 
thickness contribution to the Cauchy integral downstream of x/L = A was treated by 
assuming 6* x^/^. It is observed in Fig. 4 that the effect of including transition 

results in a smaller separation bubble and displacement thickness than that obtained 
in the completely laminar case. As a result, the edge velocity distribution in the 
transitional case more closely resembles the inviscid distribution shown in 
Fig. 4(a) since the overall interaction has been reduced in magnitude. Figure 4(b) 
shows that the skin friction has a rapid decrease at the start of transition followed 
by a rapid increase as the wall flow is energized due to turbulent mixing. This 
feature was observed in all transitional cases analyzed in the present investigation. 


Caster Experiment 

Caster (Ref. 7) performed a detailed experimental investigation on transitional 
separation bubbles which were induced on a flat plate by the pressure distribution 
which resulted from an airfoil placed near the plate. The particular case chosen 
from the Caster data set for comparison with the present analysis is denoted in 
the Caster report (Ref, 7) as Series I, number IV. The perturbation technique 
described in the analysis section was used in this calculation with the reference 
pressure, shown in Fig. 5(a), set equal to that measured by Caster with the boundary 
layer tripped from laminar to turbulent flow near the flat plate leading edge. 

The reference displacement surface was computed from a direct calculation of the 
present finite-difference boundary-layer procedure with the reference pressure 
distribution prescribed and instantaneous transition assumed at the experimental 
trip location. 

In the viscous-inviscid interaction calculation, a streamwise mesh of 81 grid 
points was distributed in the interaction region, 0.5 < s/l < 1.5, with a constant 
step size of AS/L = 0.0125, In this particular case, the best agreement of the 
computed results with the experimental data was obtained with transition assumed to 
occur instantaneously at A/L = 1.0. Figure 5 shows the predicted results for 
pressure, skin friction, displacement thickness, momentum thickness Reynolds number, 
and velocity profiles at various streamwise locations. Figure 5(a) shows excellent 
agreement between the computed pressure distribution and that obtained by Caster 
for this transitional separation bubble. The localized nature of this interaction 
is evident in Fig. 5(a) since the transitional pressure distribution merges smoothly 
with the turbulent or reference pressure distribution upstream and downstream of the 
strong interaction region. It is observed in Fig, 5(a) that the presence of the 
transitional separation bubble slightly reduces the suction peak and results in a 
pressure plateau further downstream in the separated flow region. The ’’break** in 
the pressure curve occurs at the point where the flow was assumed to undergo 
instantaneous transition with a rapid recovery in pressure taking place immediately 
downstream. Figures 5(b) and 5(c) show excellent agreement in the predicted and 
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measured separation and reattachment points and the momentum thickness Reynolds 
number at the separation point, respectively. It is observed in Figs. 5(b) and 5(c) 
that just downstream of transition, a simultaneous decrease in the displacement thick- 
ness and an increase in the momentum thickness occurs which results in a rapid decrease 
in the shape factor, H = 6 /0. In many airfoil analyses, the transitional separation 
bubble is treated in an approximate manner by assuming the shape factor decreases 
discontinuously at the laminar separation point as transition from laminar to turbulent 
flow takes place. This approximate treatment of the transition region allows the 
analysis to continue without divergence due to the presence of separation; however, 
the details of the transitional separation bubble are completely missed with this 
approximation . 

Figure 5(d) shows some of the velocity profiles computed in this case which illustrate 
the significant changes which occur in the shape of the profile over the interaction 
region. These profiles have only been plotted out to the n-location where u/u = 0.995 
in order to show the large streamwise growth which occurs in the boundary-layer 
thickness as the transition process takes place. 

Numerical tests were performed to test the sensitivity of the interaction 
analysis to the specification of the transition region- Figure 6(a) shows the two 
intermittency distributions which were used in these calculations along with the 
instantaneous transition model that was used in the results previously discussed - 
For the case designated as "small transition region" the Dhawan and Narasimha 
intermittency function given in Eq. (38) was used with Stj “ 0.9395 and ^^2 ~ 1-158; 
in the case referred to as "large transition region", the start and end of transition 
were specified by = 0,958 and 5^2 = 1.271, respectively. Figures 6(b), 6(c), 
and 6(d) show the effects of the different transition regions on the computed pressure, 
skin friction, and displacement thickness, respectively. These results show only 
small differences between the instantaneous and small transition region calculations; 
however, the computed interaction results differ significantly with the use of the 
intermittency distribution in the large transition region case. In this case, the 
onset of transition is delayed and the length scale of the transition region is 
increased as compared to the other two cases which results in a significant over- 
prediction of the extent of the separated flow as shown in Figs. 6(b), 6(c), and 6(d). 
Although the large transition region model is not applicable to this particular 
case, the computed results show the correct qualitative behavior that would occur 
if transition were delayed due to a decrease in Reynolds number. Figure 6(b) shows 
that the presence of the large separation bubble reduces the pressure suction peak; 
the pressure plateau, characteristic of separated flows, is much more pronounced 
in this case. Even though the displacement thickness is much larger in this case, 
the local interaction model is still applicable since the pressure merges smoothly 
upstream and downstream with the reference pressure distribution shown in Fig. 5(a). 
Although not presented here, results were also obtained with the imposed transition 
region moved further upstream near the laminar separation point. As expected, these 
results show a significant underprediction in the extent of the separation bubble and 
subsequent interaction. It is concluded from this study that the use of the Dhawan 
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and Narasimha intermittency function in the present interaction analysis produces 
physically correct results; however, further work is clearly required in order to 
obtain an analytical procedure which provides the onset and length of the transition 
region. 

A brief numerical study was also conducted to determine the effect of the 
streamwise mesh size on the computed results. Three calculations were made for 
the case which utilized the small transition region with the grid spacings of 
As/L = 0.025, 0.0125, and 0.00625. Figure 7(a) shows the skin friction deduced 
in these calculations in which it is observed that the results show only minor 
differences due to the change in the grid size. More sensitivity is observed in 
the computed displacement thickness shown in Fig. 7(b) in which noticeable changes 
occur downstream of the transition region for the different grid spacings which 
were used. The first order accuracy of the boundary- layer finite-difference scheme 
is apparent in Fig. 7(b) since the differences between the solutions vary in an 
approximately linear manner with the size of the grid spacing. Despite these 
differences in the displacement thickness, the computed pressure was essentially 
the same for all three grid sizes. In the future, if it is considered necessary, 
a second-order scheme could be utilized to achieve a more rapid convergence rate 
of the solution with decreasing mesh size. Alternately, as has been done in the 
remainder of the calculations, a variable grid can be used with the smallest grid 
spacings placed in the transition region to minimize the truncation errors associated 
with this high gradient region. 


Gault NACA 66^-018 Airfoil 
3 

Gault (Ref. 6) made detailed experimental measurements of transitional bubbles 
on different airfoils at various flow conditions. In the remainder of this report, 
several comparisons are made with this data using the present localized interaction 
analysis. In the first case, an analysis was conducted for the transitional bubble 
which occurs at the leading edge of a NACA 66^-018 airfoil at a chord Reynolds number 
of 1.5 X 10^ and 12 degrees angle of attack. At this angle of attack, the reference 
solution, obtained from the GRUMFOIL code (Ref. 13), was found to significantly 
overpredict the lift coefficient reported by Abbott and von Doenhoff (Ref. 35) for 
this airfoil. Correspondingly, the pressures in the strong acceleration region were 
overpredicted thereby providing an inaccurate reference solution for this case. 

These errors are probably due to the inability of the GRUMFOIL code to correctly 
model the lift decrement due to the massive trailing edge separation which occurred 
in this case at 0.75 chord on the upper surface. In order to compensate for this 
problem, the GRUMFOIL code was run at a reduced angle of attack of 11.25® which 
provided a much better match of the experimental lift coefficient and the pressure 
distribution between the stagnation point and the peak suction region as shown in 
Fig. 8(a). 

The local interaction analysis was performed with a variable streamwise mesh 
of 78 grid points distributed over the interaction region - 0.0148 < s/c < 0.31, 
where s is the arclength measured from the airfoil line of symmetry. In the present 
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case as well as the NACA 0010 (Modified) airfoil case to be discussed next, the 
arclength was measured along the actual airfoil surface Instead of the reference 
displacement surface since these two surfaces are nearly coincident for these two 
cases. The transition region was imposed from sti = 0.0258 to s ^2 " 0.0432 with 
the Dhawan and Narasimha intermit tency function. The location of st^^ was placed 
approximately midway between the separation point and the ’*break" point in the 
experimental pressure distribution; was located so that the intermittency 
function gave a value of y = 0.5 at the experimental pressure **break" point. 

With an assumed initial distribution of a linarly varying displacement thickness, 
it was found that the interaction calculation was unstable; however, this 
instability was eliminated and converged results were obtained in 50 global iteration 
cycles by using an underelaxation factor of 0.5 on the perturbation mass flow 
parameter 5 m, in each cycle of the global iteration procedure shown in Fig. 3. 

The computed results for this case are shown in Figs. 8(a) and 8(b) for the 
pressure, and skin friction and displacement thickness, respectively. Good 
agreement is obtained between the predicted results and Gault's pressure data and 
measured separation point for this case. Figure 8(a) also shows the inviscid 
airfoil pressure distribution which was obtained from the GRUMFOIL code with 
a = 11.25®. The large difference between this solution and the viscous airfoil 
solution obtained with the GRUMFOIL code shows the importance of including the 
viscous effects in the reference solution which is input to the present perturbation 
interaction analysis. 

Two additional calculations were performed to test the sensitivity of the 
interaction analysis to the upstream and downstream limits on the Cauchy integral 
given in Eq. (12). In the first calculation, the lower limit was placed at s^ = 0 
with s^ = 0.31 as was used earlier; in the second calculation, the upper limit was 
placed at s^ = 0.12 with s^ = -0.0148 as was used earlier. The results of these 
calculations resulted in no significant changes from those discussed previously, thereby 
verifying that in the present perturbation approach the numerical solution of the 
Cauchy integral can be confined to the interaction region. This point is further 
supported by Fig. 8(c) which shows the streamwise variation. of the difference 
between the interacted and reference displacement thickness. Both upstream and 
downstream of the interaction region, the increment in the displacement thickness 
is nearly constant and thus provides virtually no contribution to the Cauchy 
integral source strength in these regions. It is observed in Fig. 8(c) that upstream 
of the interaction region, in the laminar portion of the flow, the displacement 
thickness increment differs slightly from zero, this small difference being due to 
the difference between the integral laminar boundary-layer analysis in the GRUMFOIL 
code and the present finite-difference approach. 


Gault NACA 0010 (Modified) Airfoil 


Calculations were performed for this airfoil, tested experimentally by Gault 
(Ref . 6) , for an angle of attack of 8® and a chord Reynolds number of 2 and 4 x 10 
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The reference viscous airfoil solution, obtained from the GRUMFOIL code, gave good 
agreement as shown in Figs, 9(a) and 9(b) with the experimental pressure distribution 
in the strong acceleration region ahead of the transitional separation buffle. In 
this case, the turbulent separation at the trailing edge was less than 1 % chord 
thereby having no significant effect on the leading edge pressure distribution. 

In the interaction calculations, 71 grid points were distributed nonuniformly in 
the interaction region, 0 - s/c - 0.32, with the minimum grid spacing placed in 
the transition region. The intermittency functions for these calculations, which 
are shown in Fig. 9(c), were prescribed in the same manner as those used in the 
NACA 663 -OI 8 airfoil calculations. 

Comparison of the predicted pressure distributions and Gault's experimental data for 
Re^ = 4 X 10^ and 2 x 10^ are shown in Figs. 9(a) and 9(b), respectively. Although 
the present analysis shows the same overall features as the data, the detailed 
agreement between the interacted solution and Gault's data is not as good as that 
shown earlier for the NACA 66^-018 airfoil. In this case, the data has a local 
constant pressure region in the peak suction region which was not obtained in 
either the interaction solution or in the inviscid and reference solutions obtained 
from the GRUMFOIL analysis. A number of calculations were performed with different 
intermittency distributions to see if improved agreement could be obtained with the 
experimental data in the peak suction region. However, none of these assumed 
functions improved the agreement since this region is upstream of the separation 
bubble as shown by the skin friction distribution presented in Fig. 9(d). It is 
possible that this discrepancy arises due to the coarseness in the specification 
of the airfoil coordinates given by Gault for this modified NACA airfoil. 

Comparison of the pressure distributions given in Figs, 9(a) and 9(b) shows 
that the onset of the adverse pressure gradient region, which leads to the separation 
bubble, occurs further upstream in the interaction solution than that observed 
experimentally by Gault, As a result, the present laminar separation point is 
predicted upstream of that measured by Gault as shown in Fig. 9(d) for the skin 
friction distribution. In addition, the present calculation overpredicts the 
length of the separation bubble since the point of reattachment was found to occur 
downstream of that measured by Gault. Despite these differences, the present analysis 
as shown in Fig, 9(d), predicts the correct Reynolds number effect on the airfoil 
leading edge separation bubble since the separation point moves only slightly f or 
an increase in Reynolds number whereas the upstream movement of the reattachment point 
is by about the same amount as that found by Gault. Figure 9(e) shows the displacement 
thickness distributions for these two Reynolds numbers in which the upstream shift 
in the interaction region and the decreased magnitude of the displacement thickness 
are observed for the higher Reynolds number case. The pressure distributions are 
consistent with this result since the higher Reynolds number case shown in Fig. 9(a) 
has a shorter pressure plateau in the separated flow region than that for the lower 
Reynolds number case presented in Fig. 9(b). 
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CONCLUDING REMARKS 


A local inviscid-viscous interaction technique has been developed for the 
analysis of low-speed airfoil leading-edge transitional separation bubbles. In this 
analysis an inverse finite-difference boundary-layer procedure is iteratively 
combined with a Cauchy integral representation of the inviscid flow which is 
assumed to be a locally linear perturbation to a known global viscous airfoil 
solution. In this report a number of comparisons have been made with experimental 
data for transitional separation bubbles. Overall, good agreement has been obtained 
in these comparisons which indicates the validity of the present local interaction 
technique. In contrast with existing global airfoil analyses, the present localized 
interaction analysis is demonstrated to be capable of resolving the details of 
leading-edge transitional bubbles which occupy only a few percent of the airfoil 
chord. In the future this interaction analysis should be applied over a wider 
range of flows to see if the correct Reynolds number and angle of attack effects 
on airfoil transitional separation bubbles can be predicted. 

In these comparisons a forced transition model has been used in which the 
onset and length of the transition region was specified a priori. Numerical tests 
were performed with this model which indicated that the flow is somewhat sensitive 
to the specified transition region. Whereas the use of a forced transition model 
allowed the assessment of the overall interaction model presented herein, in the 
future work should be directed at the development of a natural transition model 
for the prediction of transitional separation bubbles. 
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Fig. 1 Schematic diagram of airfoil laminar-transitional separation bubble and pressure distribution 
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Fig. 4 Comparison of iaminar and transitionai separated fiow resuits 
a) Edge velocity and body shape 
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Fig. 4(b) Skin friction 
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Fig. 4(c) Displacement thickness 
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Fig. 5(b) Skin friction and dispiacement thickness 
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Fig. 5(d) Velocity profiles 
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Fig. 6 Effect of transition region on predicted resuits 
a) intermittency function 
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Fig. 8 Predicted results for NACA 663-018 airfoil 
a) Pressure 
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Fig. 9 Predicted resuits for NACA 0010 (modified) airfoii 
a) Pressure distribution for Rego^^ =4 x 10^ 
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Fig. 9(b) Pressure distribution for RSoo^c = 2 ^ 
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Fig. 9(d) Effect of Reynolds number on skin friction 
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